Gate simulation and lower bounds on the simulation time 
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Unitary operations are the building blocks of quantum programs. Our task is to design efhcient or 
optimal implementations of these unitary operations by employing the intrinsic physical resources of 
a given n-qubit system. The most common versions of this task are known as Hamiltonian simulation 
and gate simulation, where Hamiltonian simulation can be seen as an infinitesimal version of the 
general task of gate simulation. We present a Lie-theoretic approach to Hamiltonian simulation and 
gate simulation. From this, we derive lower bounds on the time complexity in the n-qubit case, 
generalizing known results to both even and odd n. To achieve this we develop a generalization 
of the so-called magic basis for two-qubits. As a corollary, we note a connection to entanglement 
measures of concurrence- type. 

PACS numbers: 03.67.Lx, 03.67.Mn 



I. INTRODUCTION 

As a starting point for the emerging field of quantum 
computation (for a review see, e.g., Refs. P, 0) Feyn- 
man Q established in 1982 a connection between simula- 
tion of quantum systems by computers on the one hand 
and computation using quantum systems on the other 
hand. Since this time many different models for compu- 
tation using quantum systems have been proposed. 

In all this models the common target is to construc- 
tively implement unitary operations by employing the in- 
trinsic physical resources. Forming the counterpart of op- 
erations on a classical computer, these unitary operations 
build the components of quantum programs. Whichever 
(reasonable) system is provided to us by an experimen- 
talist, we exploit the resources immanent to the system. 

In this paper we confine ourselves to rt-qubit systems 
where the resources are given by local unitary operations 
and by the natural time evolution specified by the Hamil- 
ton operator. This allows us to implement a given unitary 
operation by interrupting the natural time evolution with 
local unitary operations. Referring to such implementa- 
tions as programs, our objective is to study efficient or 
optimal programs. To achieve this, usually two differ- 
ent versions of this problem are considered. Hamiltonian 
simulation is the one version and denotes an infinitesi- 
mal implementation of unitary operations, i.e., the uni- 
tary operation is in a neighborhood of the identity. The 
second version is called gate simulation and describes an 
implementation of unitary operations not restricted to a 
neighborhood of the identity. We give an exact definition 
of both versions in Sec. ^1 

In this paper we address two major topics. First we 
consider two-qubit systems. Combining the manifold re- 
sults known for two-qubit systems, we can close gaps and 
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simplify the line of reasoning. In addition, by using Lie- 
theoretic methods we derive a unified approach to the 
methodology. The motivation for this extensive recon- 
sideration is twofold: First, we show that understanding 
simple cases in detail helps to generalize them to higher- 
dimensional systems. Secondly, we obtain generally ap- 
plicable (Lie-theoretic) methods which provide tools in 
the analysis of higher-dimensional systems. 

Beyond this, the two-qubit case is interesting in its 
own. We can characterize the minimal time for Hamilto- 
nian simulation directly using arguments from Lie theory. 
With the help of this characterization we explain and re- 
prove a majorization-like condition Q for the minimal 
simulation time. By employing the Weyl group of the 
corresponding Lie algebra, we are able to simplify and 
clarify the known approach of Ref. 3, especially w.r.t. 
Ref. 1^. In the case of two-qubit systems we consider 
gate simulation as well. We present a refined analysis of 
the majorization-like condition of Refs. which will 

be built explicitly on the results of Ref. Q ■ 

The second major topic deals with the case of general 
n-qubit systems. We first discuss a generalization of the 
so-called magic basis [1,0 to higher dimensional systems. 
With this information on the structure of unitary oper- 
ations we develop lower bounds on the minimal time for 
gate simulations. Our method applies to all rt-qubit sys- 
tems and generalizes a result of Ref. for even n. We 
also discuss the used techniques in connection to entan- 
glement measures, in particular, the concurrence [^lllj|. 

The whole text is written in a Lie-theoretic flavor. For 
further reference, the Lie-theoretic concepts needed in 
the text will be introduced briefly in Sec. IIIII Using this 
theory puts us in the position to formulate strong argu- 
ments in a coherent language. 

Firstly, in Sec. ^ we introduce our model. In Sec. IIIII 
we state the Lie-theoretic concepts needed in the main 
body of the text. The Hamiltonian simulation for two- 
qubits will be discussed in Sec. IIVI followed by the anal- 
ysis of gate simulation for two-qubits in Sec. The 
generalization of the magic basis for two-qubits is consid- 
ered together with lower bounds on the time complexity 
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for general n-qubit systems in Sec. I VII In Sec. IVIII we 
give a brief outline to related work, in Sec. IVIlH we con- 
tinue with a discussion of connections of our approach to 
concurrence-type entanglement measures, and in Sec. IIXI 
we close with the conclusion. In the appendix we recall 
a spectral approach to infinitesimal Hamiltonian simula- 
tion. 



II. THE MODEL 

We consider a system of n qubits, where n G N is 
finite. This system can be modeled within an n-fold ten- 
sor product (8) • • • (8) of n two-dimensional com- 
plex vector spaces C^, i.e., a tensor product of single 
qubits. The time evolution of the system is governed by 
the Schrodinger equation for the time-evolution operator, 
see, e.g., Ref. El p. 72], 

^U{t) = {-iH)Uit), 

where t denotes the time, U{t) the time evolution oper- 
ator, and H the Hamilton operator which is supposed 
to be time-independent {h — 1). Because of the irrele- 
vance of a global phase in quantum mechanics, we restrict 
ourselves to evolution operators from the special unitary 
group SU(2"). 

In addition to the possibility to let the system evolve 
according to the evolution operator U{t) = cxp{—iHt), 
in our model we allow the application of local unitary op- 
erators. A unitary operator is considered as local when it 
does not induce any interaction between different qubits, 
i.e., when it has the form of an n-fold tensor product 
Ui ■ ■ ■ Un of unitary operators Ui £ SU(2) with 
i E {l,...,n}. The time for the application of local 
unitary operations is negligible and assumed to be zero. 
Thereby we have specified the available resources which 
constitute the possibilities to control the system. 

We emphasize that, for mathematical reasons, we re- 
strict our model to consider only systems for which 
the system Hamilton operators can be represented with- 
out use of local terms. As a consequence, in the case 
of infinitesimal Hamiltonian simulation we can simulate 
Hamilton operators exactly only if the Hamilton oper- 
ator can be represented without use of local terms. To 
remove the local terms of a system Hamilton operator one 
usually employs some approximations, see, e.g., Ref. 0, 
p. 3] or Ref. [Ij, p. 288]. But we refrain from consid- 
ering such approximations. In addition, avoiding local 
terms in Hamilton operators seems to release us from 
some problems with infinite programs, i.e., an infinite 
number m = oo of steps (see below). In Ref. 0] it is 
analyzed under which conditions we need infinite pro- 
grams for time optimal control in the more general set- 
ting which includes local terms in the system Hamilton 
operator. Nevertheless, we consider for technical reasons 
all types of programs, even infinite ones. The available 
resources will be utilized below in three ways. 



First, we consider the simulation of a unitary gate, i.e., 
a unitary operator. This notion of simulation means that 
the system is able to implement a given unitary gate by 
interrupting the natural time evolution with local unitary 
operations. Dcfinition^states this in a more formal wa; 
The term "gate simulation" was introduced in Ref. 
p. 3]. 

Definition 1 (Gate simulation). An n-qubit system 
with Hamilton operator H and local unitary operators 
available simulates a unitary gate U in time t if there 
exists local unitary operators Uo and Uj as well as times 
tj > with t — j G r^, 1 < j < and m G 

NU {0,cx)}, so that 



U 



Uo. 



Remark. Due to the non-commutativity of the unitary 
group we restrict the symbol Y[ for elements Vj from the 
unitary group to the following meaning: 



n for / > e, 

id for / < e. 



where e, / G Z, and the identity element of the unitary 
group is denoted by id. In the case of / infinite the 
symbol Y[ represents an element from the closure of con- 
vergent sequences. 

Secondly, we introduce a particular concept of in- 
finitesimal simulation of a unitary gate. A given uni- 
tary gate U is treated as a point of a one-parameter 
group exp(— iiJ't'). The infinitesimal simulation of the 
Hamilton operator H' denotes that the system simulates 
the corresponding one-parameter group for infinitesimal 
times t' , i.e., the derivatives of the one-parameter group 
and of the simulation coincide for infinitesimal times. 
We emphasize that in Def. |21 the notion of infinitesimal 
Hamiltonian simulation is defined independently of the 
unitary gate U. 

Definition 2 (Infinitesimal Hamiltonian simula- 
tion) . An n-qubit system with Hamilton operator H and 
local unitary operators available simulates an Hamilton 
operator H' infinitesimally in time t if there exists local 
unitary operators Uq and Uj as well as times tj > where 
* = Ejli j € N, 1 < j < m, and m € N U {0, oo} 
so that njLo Uj is equal to the identity of the unitary 
group and the following equation holds: 



lim 

t'— 
t'>o 




Y[ Ujexp{-it'Htj) 



Un 



(1) 
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Remark. The condition 11^=0 ~ ^'^ ensures that our 
program specified by Uq, the UjS, and the tjS operates 
nearby the identity for t' . 

We note that in Refs. [lllil|ll[l3| the notion of in- 
finitesimal Hamiltonian simulation was extended by the 
so-called first order approximation to unitary operators. 
Similar ideas were used in Ref. Q where it was pro- 
posed to follow the evolution of the system exactly by the 
Hamiltonian simulation. Concurrently, it was remarked 
in Ref. 4] that to follow the evolution of system exactly is 
only possible infinitesimally, as the control is not continu- 
ous. In this text we do not consider such approximations. 

Though Def. |2] presents the essential meaning of in- 
finitesimal Hamiltonian simulation, it appears to be very 
unpractical. Thus we present an equivalent condition, 
which is usually formulated as definition [Tsl IT^ IT^ 

Ha]. 

Lemma 1. An n-qubit system with Hamilton operator H 
and local unitary operators available simulates the Hamil- 
ton operator H' infinitesimally in time t if and only if 
there exists local unitary operators Vj as well as times 
tj > with t — "^j^itj, j (z N; 1 < j < m, and 
TO G N U {0, oo} so that the following equation holds: 



Remark. Actually, we do not use Def. plater in the text. 
But we state this definition to highlight that Def. |21 is 
independent of some unitary operator U and does not 
incorporate different decompositions of U which could 
lead to different Hamilton operators H' . The existence 
of different Hamilton operators H' will be employed in 
Sec. El below. 

Before we proceed, we discuss our model. Entangle- 
ment describes important non-local properties of states 
and gates. Because entanglement is invariant under lo- 
cal unitary operations it seems reasonable to neglect 
the time needed to implement local unitary gates for the 
implementation of general unitary gates. This is sup- 
ported by the fact that two-qubit gates are considered as 
significantly more difficult to implement than one-qubit 
gates 20] . Additionally, in Nuclear Magnetic Resonance 
(NMR) the application of local unitary operations in zero 
time is expressed by the notion of the "fast control limit" 
which is conventionally considered as a good approxima- 
tion. This is reasonable because local and non-local gates 
operate on different time scales |^ |^ . 

III. LIE-THEORETIC PREPARATIONS 



3 = 1 



(2) 



Proof The "only if'-case: The l.h.s. of Eq. (P) equals 
—iH' and the r.h.s. of Eq. equals 



lim 

t'>0 



where 




Y[exp[-it'{WjHWj-^)tj] ^Wo 

i=i 



(3) 



U„ 



for j 



Wj+iUj for < j < TO. 



We differentiate, compute the limit, and equate the result 
of Eq. lO with ~iH'. After this we use Vj := Wj^^ and 
Wo = njLo = id to obtain Eq. @. 

The "if'-case: After insertion of Eq. (0) in Eq. ^ we 
obtain the "if'-case. □ 

Thirdly, we introduce in Def 13 the notion of infinitesi- 
mal gate simulation which depends explicitly on the given 
unitary gate U. 

Definition 3 (Infinitesimal gate simulation). An n- 

qubit system with Hamilton operator H and local unitary 
operators available simulates a unitary gate U infinites- 
imally in time t if there exists an Hamilton operator H' 
and local unitary gates Ui and U2 such that the equation 
U = Ui exp{—iH')U2 holds and the system simulates the 
Hamilton operator H' infinitesimally in time t. 



In this section we recall Lie-theoretic notions and 
methods which will be employed throughout the paper. 
This reflects the intimate connection of the considered 
problems to Lie theory and it makes the text more read- 
able and self-contained. For convenience, this section can 
also be regarded as a reference section. We remark that 
our presentation in this section was partly inspired by 
Refs. 



A. Basic concepts 

For our purposes we can consider Lie groups as linear 
matrix groups, i.e., as closed subgroups of the general 
linear group. The tangent space to the Lie group G at the 
identity is isomorphic to the Lie algebra g corresponding 
to G. A Lie algebra, which is in particular a vector space, 
comes with a bilinear and skew-symmetric multiplication 
operation called the Lie bracket [ , ] . The Lie algebra g 
is closed under the Lie bracket and the Jacobi identity 

[[31, 52], 53] + [[33, 51], 52] + [[32, 53], 51] = 

holds for all elements 31,52,53 G 0- We emphasize that 
we use only real or complex Lie algebras which are finite- 
dimensional. For general reference on Lie groups and Lie 
algebrasplease consult Refs. jH iilil [il 111 HI IH 
IS EE IS In the following let G denote a Lie group 
and 2 its associated Lie algebra. 

The map adg{g) from the Lie algebra g to itself is 
defined by /i 1-^ [g,h], where h,g S g. With this 
notation, the adjoint representation of the Lie algebra 
g in itself is given by 5 1-^ a,dg{g). Let adg(f)) de- 
note the set {adg(/i)|/i G f)} for some subspace t) of 
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the Lie algebra g. Now, we can introduce a symmet- 
ric bilinear form on the Lie algebra g: the Killing form 
Bg{g,h) := Trg(adg((7) o adg(/i)). If the Killing form is 
non-degenerate it can be thought as an inner product on 
the Lie algebra, although it does in general not fulfill the 
axiom of positivity. 

Definitio n 4 (Orthogonal symmetric Lie algebra, 
see Ref. [25', p. 213] and Ref. [30, pp. 225 226 
and 246]). A pair (g, 6) is an orthogonal symmetric Lie 
algebra if 

(i) g is a real Lie algebra, 

(ii) 9 is an involutive automorphism of g, 

(iii) and the connected Lie group of linear transforma- 
tions of g generated by adg(6) is compact, where 6 
is the set of fixed points of 6* in g. 

Remark. An automorphism of a Lie algebra g respects 
the Lie bracket, i.e., for all g and /i in g we have 0{[g, h]) — 
[9(g), 9{h)]. An involutive automorphism is in addition 
self-inverse. Let i and p be the eigenspaces of 61 in g 
for the +1 and —1 eigenvalue respectively. Consider the 
canonical decomposition g ~ i + p. Condition (ii) of 
Def . ^ is equivalent to 

C i, [t,p] C p, [p,p] C i (see [13, p. 226-227]). (4) 

In Ref. 5] this decomposition was called Cartan decom- 
position. If Eq. Q holds we can define 9 by 

9{k) = k for &\\k el and 9{p) = -p for all pep. (5) 

In addition when g is the Lie algebra of a compact group 
G then Condition (iii) of Def. 0] is always true. 

For further reference, we assume that g is semisimple, 
i.e., that the Killing form of g is non-degenerate, and 
that (g, 6*) is an orthogonal symmetric Lie algebra. As 
an example, the Lie algebra su(2"), which corresponds to 
the Lie group SU(2"), is semisimple. We fix a canonical 
decomposition g = 6 -I- p satisfying Eq. Q and a maximal 
Abelian subalgebra a contained in p. Let K — exp(t) and 
A = exp(a) denote the subgroups of G generated by 6 
and o respectively. After this preparation we obtain a 
decomposition G = K AK of the Lie group G. 

Fac t 1 (K A K decomposition of the Lie group G 
[25', Ch. V, Thm. 6.7]). With the notation as given 
above, the Lie group G corresponding to g can be decom- 
posed as 

G = KAK. 

Similar to the adjoint representation adg of a Lie al- 
gebra g in itself, we can define the adjoint representa- 
tion Adg of a Lie group G in its Lie algebra g. For 
an element G G G we introduce 4>g{G) as the map 
H G~^HG with the signature G ^ G. The map 
Adg(G') has the signature g ^ g and is defined as the 



differential of (/>g(G). For matrix representations we can 
write Adg(G') as the map g ^ G~^gG. We use the short- 
cut Adg(K) := U^gj^ Adg(iir) and get the relation be- 
tween the subspace p and its Abelian subalgebra a: 



Fact 2 ([25, Ch. V, Lemma 6.3 (iii)]). The following 
equation holds: 

p = (Adg(K))(a). 



B. The Weyl group and infinitesimal convexity 

We use the notation 
CK(a) := {K e K| (Adg(X))(a) = a for all a e a} 
and 

iVK(a) := {K e K| (Adg(if))(a) C a} 

respectively for the centralizer G]^{a) and the normalizer 
iVK(a) of a in K. 



Definition 5 (Weyl group, see |25|, p. 284] or |27 
p. 381]). The Weyl group corresponding to a is the 
factor group A'^K(a)/CK(o). We denote this group by 
>V(G, A), where A exp(o). 

The Weyl group >V(G, A) is finite (see Fact 13 below) . 
In order to compute the Weyl group, we introduce the 
concept of restricted roots. 

Definition 6 (Restricted root, cf. |,27fc p. 370]). 

Let A be a linear function on a. The linear subspace g^ 
is given by 

g'*' = {g e g| [a, g] = \{a)g for all a e a}. 

The linear function A is called a restricted root of g w.r.t. 
o if g^ ^ {0} and A is not identically zero on a. Let Ao 
denote the set of restricted roots of g w.r.t. a. 

Remark. In Ref. [23, p. 370] the restricted roots are de- 
fined w.r.t. ia. But the concept of restricted roots can 
also be defined w.r.t. a. 

Due to the fact that g is semisimple we deduce that 
the Killing form restricted to a x o is non-degenerate. 
With this in mind, a restricted root A is equal to the 
map a 1— > B^{a\,a), where ax e a is uniquely deter- 
mined. We extend the Killing form to restricted roots by 
Bg{X,^) := Bg{a\,a^i). For every A € An the refiection 
sx{p) of a restricted root fi e A^ w.r.t. the hyperplane 
{a e a\ A(a) = 0} is given by 



sa(m) 



Sg(A,A)- 



A. 



Following Ref. (25i . p. 286] the refiection s\ can be ex- 
tended to elements of a. For a e a the reflection s\ (a) of 
a G a in the hyperplane {a e a\ A(a) = 0} is given by 



s\ia) 



Bg{a,ax) 
Bg{ax,ax) 



(6) 
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With this preparation we get a possibihty to compute 
the Weyl group corresponding to a: 

Fact 3 ( j27L p. 383]). The Weyl group corresponding 
to a is finite and is generated by the reflections s\, where 
A e Aa. 

Recall that W(G,A) is a subset of K and operates on 
a by Adg{K), where K <eK. 

Definition 7 (Weyl orbit, see, e.g., fsd, p. 422]). 

The Weyl orbit W{a) of a G a is defined as the set 
{{AdgiW)){a)\W e W(G,A)}. By appealing to FactH 
the Weyl orbit yV(p), for p € p, is defined as yV(p) := 
yV{a), where a £ {Adg(K)){p) n a. 

To understand that the definition of W(p) for p e p is 
independent of a and thus well defined we now charac- 
terize the Weyl orbits in more detail. 

Fact 4 ([27, Lemma 7.38]). Let {Adg{K)){a) = a', 
where a, a' £ a and K £ K. Then there exists an element 
K' e iVK(a) such that {Adg{K')){a) = a'. 

By Fact 0] two elements a and a' from the Weyl orbit 
yV(p) of p € p are conjugated by an element of the Weyl 
group, which proves that the definition (Def. [7|l of yV(p) 
is independent of a. This shows in addition that the Weyl 
orbit W{p) is equal to an (Adg(K))(p). Let us denote the 
convex hull of the Weyl orbit W{p) by t[p). We state now 
the infinitesimal version of Kostant's convexity theorem. 

Fact 5 (Kostant's convexity theorem (infinitesi- 
mal version), see [36, Thm. 8.2] or [37, Thm. 1]). 

Let r he the orthogonal projection of p on a w.r.t. to the 
Killing form. For every p £ p one obtains 

r((Ad,(K))(p)) = c{p). 

Remark. The essential meaning of the infinitesimal ver- 
sion of Kostant's convexity theorem (FactlSjl is that the 
projection of (Adg(K))(p) to a w.r.t. the Killing form is a 
convex set and its extreme points are given by the Weyl 
orbit Wij)). 

In order to characterize the Weyl orbits in more de- 
tail we introduce additional concepts. The subspace a 
can be divided into connected components called Weyl 
chambers. 

Definition 8 (Weyl chamber (H, p. 287]). Let A 

be any restricted root of q w.r.t. a. The hyperplanes 
{a E o| A(a) — 0} divide o into finitely many connected 
components excluding their boundary hyperplanes. Such 
a connected component is called Weyl chamber. The clo- 
sure of a Weyl chamber including their boundary hyper- 
planes is called closed Weyl chamber. 

The Weyl group acts on the Weyl chambers: 

Fact 6 ([^2?, Thm. 2.12, Ch. VII]). The Weyl group 
permutes the Weyl chambers. 



We choose some arbitrary, but fixed, order on the re- 
stricted roots of 2 w.r.t. a. Therefore, the restricted roots 
can be divided into positive and negative (restricted) 
roots, where positive and negative is defined regarding 
to the chosen order. A restricted root is called funda- 
mental if it is positive and not a sum of two positive 
(restricted) roots [2^ p. 59]. Let C be the 

set of fundamental (restricted) roots. Since the Killing 
form restricted to a x o is non-degenerate we can de- 
fine as before for every (restricted) root A the element 
flA e a so that Bg{a\,a) — A(a) for all a & a. The set 
{a £ a\ Bg{aak,a) > for all at} is a Weyl chamber and 
it is called the fundamental Weyl chamber 0, p. 61]. 

Fact 7 (adapted from HE Prop. I, Sec. 2.11]). Let 

A and fi be restricted roots corresponding to elements ax 
and of the closed fundamental Weyl chamber, respec- 
tively. The element lies in the convex hull of the Weyl 
orbit W{a\) of a\ if and only if X{a) > fJ.{a) for all el- 
ements a of the fundamental Weyl chamber. The condi- 
tion A(a) > /i(a) is equivalent to Bg{a\,a) > _Bg(a^,a). 

C. The two-qubit case 

We now treat the case G = SU(4). To be more con- 
crete we introduce a matrix representation for the real 
semisimple Lie algebra su(4) which corresponds to G. 
Let 

<J. := (I J) , ay := , and a, := (^J 

be the the Pauli matrices and set 



to be the identity matrix. We identify cf^ = ai, (Jy = (72, 
and az = cr^ and use the following definitions 







gl CTl , ^2 




D a2,X3 


2-0 « 


' crs, 


Xi 


i 


S) cro,X5 


i 


5 (To, Xq :— 




' -0, 


Xj 


i 


g) (Tl , 


i 

'■— — CTt $ 

2 


D (72, ATg 




' 


Xio : 


i 


'^(72,Xii 


i 


g CF3,Xi2 := 


i 


(Xl (Tl 


Xi3 : 


i 

:^-a2^ 


5 ^3, Xii 


i 


g) CTi, := 


i 

- —a-i 
2 ^ 


(8l (72 



The standard (or defining) representation of 5u(4) is 



:=5u(4) = span^lXi, . . . , ATig}, 

where span^ denotes the real span. Let 

e := spanjj,{Xi, . . . ,^6}, 
p := spanjj{X7, . . .,Xi5}, 
and a :— span^{Xj, . . . , Xg}. 
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With this notation one can easily check that t and p fulfill 
the commutator relations in Eq. Q). Since the group 
SU(4) is compact, the pair (g, 6) defines an orthogonal 
symmetric Lie algebra, where is given by Eq. jSJl- The 
subspace o forms a maximal Abelian subalgebra in p. 
The set of restricted roots w.r.t. a can be computed as 
the eigenvalues of a.dg{ciXi + C2Xs + c^Xg): 

{±i(c2 - C3), ±i(c2 + C3), ±i(ci - C3), 

± i(ci + C3), ±i(ci + C2), ±i(ci - C2)}. (7) 

We use Eq. © to obtain a generating set for the Weyl 
group (corresponding to a) as a set of matrices 





1 
-1 



which operate on the vectors 

Oj ^Xj, llj ^Xs, and f ) ^ Xg 



, (8) 



(9) 



With the notation of Eq. the Killing form restricted 
to a X a is given by 



BJa,b)\ 

fl ^ ' '\axa 



Now we give the elements ax €z a corresponding to the 
restricted roots A in Eq. (0, i.e. elements oa G a such 
that Bg{a\, a) ~ A (a) for all a e a: 





We have used the basis of Eq. @ to represent the ele- 
ments ax- 

In order to present our results in the context of Ref. jj] 
we choose an order on the (restricted) roots, such that 
the roots of Eq. iQ which have a plus sign constitute the 
positive ones. With this convention for an element d — 
diXr + d2Xs + d^Xg of the fundamental Weyl chamber 
we get the set of equations 

{d2 -d3>0,d2+d3> 0, di - d3 > 0, 

di+d3> 0, di + da > 0, di - da > 0} , (10) 

where we have identified (M, >) with (iM, >) by defining 
iri > ir2 <4> ri > r2 for all ri, ra € M. 



IV. INFINITESIMAL HAMILTONIAN 
SIMULATION FOR TWO QUBITS 

A. Lie-theoretic explanation 

Following Def. |21we consider now infinitesimal Hamil- 
tonian simulation for two qubits. We emphasize that in 
the two-qubit case local unitary operations correspond to 
elements of K = exp(6). We use the notation of Sec. lIIII 
especially that of Sec. IIII CI Since we restrict ourselves 
to Hamilton operators without local terms (see Sec. Ill|. 
we have for all non-local Hamilton operators H and H' 
that iH £ p and iH' g p, where p is the subvector space 
of the Lie algebra g introduced in Sec. IIIII Thus, we can 
use Fact|21to write every non-local Hamilton operator H' 
as H' — (Adg((L')~^))(a'), where a' is an element of 
and L' is a local unitary operator. 

Theorem 1. Assume that H and H' are non-local 
Hamilton operators acting on a two-qubit system. Let a' 
be an element of a, where a' ~ {Adg{L')){H') for some 
local unitary operator L' . 

A two-qubit system with Hamilton operator H and local 
unitary operators available is able to simulate the Hamil- 
ton operator H' in time t if and only if the Hamilton op- 
erator (a'/t) lies in the convex closure of the Weyl orbit 
yV(iJ) of H. The condition is independent of the choice 
of a'. 

Remark. Actually, Thm. ^ is the infinitesimal version of 
Fact IHl (see below and Ref. ) . In order to clarify the 
connection of Thm. [21 to the work of Ref. |^ , we give 
here a proof of this infinitesimal version using arguments 
of Refs. HH. 



Proof. Assume that 



tj2qnAd,iKnm) 



(11) 



is a simulation of H' in time t, where K^' are local unitary 
operators, H,H' G p, q'/ > 0, and J2Z\ I'l = 1- 

Due to Fact [21 there exists a and a' in a, where 
a = iAdg{L)){H) and a' = (Adg(L'))(ff') for some lo- 
cal unitary operators L and L'. We remark that local 
unitary operations cost no time. Thus, the existence of 
the simulation in Eq. Hll|l is equivalent to the existence 
of a simulation i 9i(^dg(i4r|))(a) of a' by a in time 
t, where K'^ are some local unitary operators, > 0, 
and J2i^i 1i — 1- r and F' denote the orthogonal 
projections (w.r.t. the Killing form) of p on a and 
respectively. We can write the simulation as 

m2 

tJ2ql [r((Ad,(ifD)(«)) +r'((Ad,(ifD)(«))] =«'• 

This is equivalent to 

m 

t5^g,r((Ad,(if,))(a)) =a', (12) 
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where Ki are some local unitary operators, qi > 0, and 
Sl^i 9i ~ ^- "^^^ ^^^^ equivalence follows in the back- 
ward direction by employing FactElto rewrite the projec- 
tion r as a convex combination and in the forward direc- 
tion by the fact that the term J2TJi 9*r'((AdB(i^;))(a)) 
has to be zero. 

By the remark following Fact|Slwe know that the pro- 
jection of (Adg(K))(p) to a w.r.t. the Killing form is a 
convex set. Thus, we can write Eq. 112|l as 

r((Ad,(if'))(«)) = (a'/t) 

for some local unitary operator K' . With Fact [3 we get 
that {a'/t) lies in the convex closure of the Weyl orbit 
W(a) of a. Since a — {Adg{L)){H), we can replace a by 
H in the preceding sentence. This proves the theorem 
except for the independence of the choice of a'. 

Assume that we replace a' by a" g a, where a" = 
{Adg{L")){H') for some local unitary operator L" . Due 
to Fact^ there exist an element W G W{G,A) so that 
a" ~ (Adg(VF~^))(a'). Since operating with an ele- 
ment of the Weyl group leaves the Weyl orbit yV(7?) un- 
changed, the Weyl orbit is equal to {Adg{W)){W{H)). 
It is obvious that the convex closure of the Weyl orbit 
yy(iJ) is left unchanged as well. Hence, the element a' 
is in the convex closure of the Weyl orbit W{H) if and 
only if a" is. □ 

For a' e a it was also proven in Ref. that the set of 
Hamilton operators {a'/t) which can be simulated in time 
one is convex. We emphasize that the extreme points of 
this set are given by the Weyl orbit >V(iJ), which can be 
computed by means of Eq. ©. In Ref. U] the extreme 
points were given and their extremality was proven by 
another method. As in Ref. we state now a version of 
Thm. n which gives a condition for infinitesimal Hamil- 
tonian simulation in the two-qubit case that is easier to 
check. 

Theorem 2 (f3, p. 11]). Assume that H and H' 
are non-local Hamilton operators acting on a two-qubit 
system. Let a and a' be elements of the closed fun- 
damental Weyl chamber, where a — aiXj + ci2Xg, -\- 
asXg = {Adg{L)){H) and a' = a'^X-j + al^X^ + a^jATg = 
(Adg(L'))(iJ') for some local unitary operators L and L' . 

A two-qubit system with Hamilton operator H and local 
unitary operators available is able to simulate the Hamil- 
ton operator H' in time t iff the following equations hold: 

ai > a'l/t, (13a) 
ai -\- 02 -\- as > {a'l -I- aj, -f a's)/t, (13b) 
0-1 + 02 — as > {a'l -\- a'2 — a's)/t. (13c) 

Remark. We force a and a' to be (almost) unique ele- 
ments of a by choosing them to be elements of the closed 
fundamental Weyl chamber. If a or a' lies on the bound- 
ary of the fundamental Weyl chamber, they are elements 



of the closed fundamental Weyl chamber, but not ele- 
ments of the fundamental Weyl chamber. Only in this 
case there remains a non-uniqueness, and the considered 
element can possibly chosen to lie on different boundary 
hyperplanes of the closed fundamental Weyl chamber. 

Proof. As the Weyl group permutes the Weyl chambers 
(see Fact O we can choose a and a' to be elements of 
the closed fundamental Weyl chamber. We recall from 
Eq. ifTUI) that an element d — diX-j -\- d2Xs -\- d^Xg 
lies in the fundamental Weyl chamber iff o?2 — ^3 > 0, 
d2 + ds > 0, di - ds > 0, di -t- ds > 0, di -I- ^2 > 0, 
and di — c?2 > holds. Applying Thm. ^ and Fact [7| we 
get that oidi -|- 02^2 + 03^3 > (a'ldi + 02*^2 + 03^3)/^ 
holds for all elements d — diXj -\- d2Xs -{- d^Xg of the 
fundamental Weyl chamber. Eliminating the quantifiers 
in the previous condition, e.g ., using the computer alge- 
bra system QEPCAD (38L l39l |. we obtain the conditions 
of Eq. inj. □ 

B. Majorization 

In this subsection we introduce some concepts from the 
theory of majorization which will be employed later. Our 
presentation is succinct and we refer to Refs. 1^ 
113,13143 for a more detailed treatment of this topic. 

For an element x = {xi, . . . , Xk)'^ of R*^ we denote by 
x^ — {x{, . . . , x'l)'^ a permutation of x so that x'l > if 
i < j, where I < i, j < k. 

Definition 9 (Majorization | 45l p. 28]). A vector 
a; S R*^ is majorized by a vector ?/ G R*^ if 

'^^l <Y1 all 1 < / < fc 

i=l 1=1 

and 

k k 
i=l i=l 

The notation x ^ y means that x is majorized by y. 

We recall the notion of s-majorization introduced in 
Ref. 3. For an element x — {xi, X2, xj)"^ of R^ we intro- 
duce the vector x = {\xi\, \x2\, |a;3|)'^, and we define the 
s-ordered version x^" ofx by setting xl'' :— x{, X2' '■= £2, 
and := sgn(a::ia;2a;3)a;3. The signum of X1X2X3 is de- 
noted by sgn(a;ia;2a;3). 

Definition 10 ([4, p. 11]). The vector x £ R'"^ is s- 
majorized by ?/ € M'^ if 

xi" < yi", 
+ + xl; < y[' + + y^" , 
xi" + x^2' - x^" < yi" + y^2' - yl" . 
The notation x <a y means that x is s-majorized by y. 
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We emphasize that a vector representing an element 
from the Lie subalgebra a is s-ordered if and only if it 
lies in the closed fundamental Weyl chamber, as given 
in Eq. H1U|I . except that for the closure the relation < 
has to be replaced by the relation <. This gives a geo- 
metric motivation for the s-ordered vectors. In addition, 
the necessary and sufficient conditions for Hamiltonian 
simulation in Eq. H13|l are equivalent to the definition of 
s-majorization. 

Corollary 1 ([4, p. 11]). Assume that H and H' 
are non-local Hamilton operators acting on a two-quhit 
system. Let a and a' he elements of a, where a — 
aiXr + a2Xs + 03X9 = (Adg(L))(i?) and a' = a[X7 + 
02X8-1-03X9 = (Adg(L'))(_ff') for some local unitary op- 
erators L and L' . We use the notation a = (01,02,03)^ 
and a' = (a'j^, a'2, 03)"^. 

A two-quhit system with Hamilton operator H and local 
unitary operators available is able to simulate the Hamil- 
ton operator H' in time t if and only if the following 
equation holds: 

a' td. 

There is a similar condition on infinitesimal Hamilto- 
nian simulation which is given in terms of a majorization 
condition on the spectra of the considered Hamilton oper- 
ators. This result (see appendix and Ref. ^id, pp. 9-10]) 
should be compared to the s-majorization condition in 
Cor.Hl 

V. GATE SIMULATION FOR TWO QUBITS 

As in Def. ^we consider now gate simulation which is 
a global version of infinitesimal Hamiltonian simulation. 
We recall a theorem of Khaneja et al. |^]. 

Fact 8 ([5, Thm. 10]). Assume that H is a non-local 
Hamilton operator acting on a two-quhit system. 

A two-quhit system with Hamilton operator H and local 
unitary operators available is able to simulate the unitary 
gate U in time t if and only if the unitary gate U can be 
decomposed as 

(7 = iiexp(W)i2, (14) 

where Li and L2 are local unitary operators and W is an 
element which lies in the convex hull of the Weyl orbit 
W{H) ofH. 

Remark. An equivalent version of Eq. 114|) is 

L-if/L-i = exp(tW^). (15) 

This means that U can be simulated in time t if and only 
if there exists a unitary gate U' which is locally equiva- 
lent to U and which can be expressed as U' = exp{tW). 
But there exists a restriction on the elements Li and 
L2. As exp{tW) is an element of A = exp(a), we have 



that L^^UL2^ has to be an element of A, too. There 
exists different unitary operations U' which satisfy this 
restriction. The appearance of different unitary opera- 
tions U' is a consequence of the non-uniqueness of the 
K A K decomposition of Fact ^ which will be analyzed in 
detail below. We emphasize that it may be impossible 
to express U as U = exp{tW) with the same (or shorter) 
time t as in Eq. (|15|) . 

We present now the results on gate simulation in simi- 
lar fashion as done in Section If V Iji for Hamiltonian sim- 
ulation. Due to the remark following Fact |S| a local uni- 
tary operation U can be simulated in time t if and only if 
a local unitary operation U' which is locally equivalent to 
U can be expressed as U' — exp{tW), where W denotes 
an element of the Weyl orbit of the system Hamiltonian. 
In the sequel, let Ki, for i g {1, . . . ,8}, be suitable el- 
ements from the set of local unitary gates K = exp(t). 
In addition we denote by A and A' some appropriate 
elements of A = exp(a). In view of Fact ^ we can de- 
compose the unitary gate U and locally equivalent gates 
U' asU = K1AK2 and U' = K^A'Ki respectively. To 
characterize all unitary gates U' which are locally equiv- 
alent to U it is necessary and sufficient to characterize 
all A' satisfying K^A'Kq = A. Thus we have to identify 
all A' which can be written as A' = [K^^ AK'j)Kg,. This 
is done in the following lemma. 

Lemma 2. For a fixed A e A and an arbitrary element 
of the form A' = {K-^AK)K' e A with K,K' &K we 
can choose K from the Weyl group and K' from the set 

Kn A. 

Related to 9 from the definition (Def . 0J of an orthogo- 
nal symmetric Lie algebra (g, 9), there exists a global ver- 
sion 9 operating on the Lie group G, see, e.g., Ref. [s^, 
Thm. 2.3 of Chap. IV] or Ref. [H Thm. 6.31.]. We de- 
fine e by Q{K") = K" for K" e K and e(P) = for 
P £ P — exp(p). We employ the mapping ()*: G ^ G 
given by G 1-^ G* := 6(G^^). We use here the sym- 
bol 0* in order to avoid confusion with the symbol 
()*, which denotes complex conjugation. We have that 
(G'iG2)* = G*G* for Gi,G2 e G, P* = P for P e P, 
and {K"y = {K")-'^ for K" e K, see Ref. [M p. 81]. 
We introduce the map (j>: G/K P which is defined as 
GK i-> (j)iGK) ■.= (GK){GK)* = GG*. This map cj) was 
studied in Ref. [lip. 81-82] and in Ref. [H Proof of 
Thm. 6.31.]. Ref. jSJ] proves that (j) induces an isomor- 
phism of G/K onto P. 

Proof of Lemma\^ We employ the map (f) and obtain 
the equations (j){{K-^ AK)K') = R-^A^K and (^(A') = 
{A')'^. Since A' = {K~^AK)K' is given in the condition 
of Lemma 13 we obtain K'^A^K = {A')'^. And due to 
Facial we can choose K as an element of the Weyl group. 
Thus, K^^AK e A which proves that G K n A. □ 

We still need to characterize the elements of K n A. 
This will be done now. 
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Lemma 3. The elements of set K n A are given by 
exp{ziTTXi + Z2TrX^ + z^ttXq), where Zj € Z for j G 
{1,2,3} and Xy, Xg, Xg as defined on page\^ 

Proof. First, we show that the elements exp(zi7rX7 + 
Z2TtXs + z^nXg) constitute a subset of K n A. Since 
e'x.'p{ziTrXT + Z2TrXs + zsttXq) for Zj G Z are by definition 
elements of A and A is an Abelian group, we obtain that 

exp(zi7rX7 + Z2nXs + zsTrXg) 
— exjp{ziTT Xt) exp(z27rAr8) cxp(z37rXg) 
= (icri ig) aiy^{ia2 «> cr2)^H«c^3 «> (JsY^ . 

This proves that the elements constitute a subset of KnA. 
Secondly, we show that K n A is a subset of the set given 
by the elements exp(zi7rA"7 + Z2TtXs + z^ttXq). We make 
the ansatz exp{ajXj+ag,Xg,+agXQ) = cxp(aiXi+a2X2+ 
asX^ + aiXi+a^sX^ + aeXe), where € M, Xi were given 
on p. |31 and i G {!,..., 9}. By direct computations one 
gets for a7, os, and ag the conditions 

(ai — as — o,q)/tt G Z, 
{ar + as — ag)/TT e Z, 
{ar + as + ag)/TT e Z, 
(a7 — as + ag)/n G Z. 

This implies ai/ir e Z for « e {7, 8, 9}. □ 

Now, we state the majorization-like equivalent of 
FactEl 

Corollary 2 (see Ref. [6, Lemma] or Ref. 
Result 1]). Assume that H is a non-local Hamil- 
ton operator acting on a two-qubit system and that we 
intend to simulate the unitary operation U. Let a 
and a' be elements of a, where a — aiXi + a2Xs + 
03^9 = {Ad^{Ki)){H), a' = a'iX7 + ^Xg + a'^^Xg, 
and U = K2 exp(a')i\r3 for some local unitary operations 
Ki, K2, K3 G K. We use the notation a — (01,02,03)"'" 
and a' — {a[, a'2, 03)"^- 

A two-qubit system with Hamilton operator H and local 
unitary operators available is able to simulate the unitary 
operation U in time t if and only if the following equation 
holds for at least one choice of z — (21,22,-^3)"^ G 

a' + TTZ -<s ta. 

Proof. By Fact |21 and Fact we can choose a and a' 
respectively as given. Applying the remark following 
Fact |S1 it is necessary and sufficient to consider some 
unitary gates U' which are locally equivalent to U. By 
use of Fact ^ these locally equivalent gates U' can be 
represented as U' — K[A'K2, where A' is an element 
of A and K[, K2 are local unitary gates. The differ- 
ent possibilities for A' in this decomposition are given 
by Lemma 12 as A' = exp[{Adg{K)){a')] K' , where K 
is an element of the Weyl group, K' G KnA, and 
K' = exp(fc'). With the characterization of K n A 
from Lemma |21 we deduce that K n A is left invariant 



by operations of the Weyl group. Since A is Abelian 
and K n A is left invariant by operations of the Weyl 
group, we can write A' as A' — exp [{Adg{K)){a') + k'] — 
exp[{Adg{K))(a' + k")], where K" = exp(fc") for some 
element K" G K A. By Fact we obtain that A' = 
exp [{Ad^{K)){a' + k")] = exp(W), where W lies in the 
convex hull of the Weyl orbit W(a) of a. When we con- 
sider the equation exp [{Adg{K)){a' -\- k")] = exp{tW) in 
a basis where both (Adg(if ))(a' -I- k") and tW are diago- 
nal then we obtain by the periodicity of the exponential 
function that {Adg{K)){a' k") -\- M = tW, where M = 
diag(27riAi, 27riA2, 27rjA3, 27riA4) and Ai,A2,A3,A4 G Z. 
Since {Adg{K)){a' -\- k") and tW are elements of a it 
follows that M G a. We can write M as M = 2nziX'! -\- 
2'KZ2Xs -f 27123X9 = 2ki where zi = (Ai -I- A2) G Z, 
Z2 = (A1-I-A3) e Z, Z3 = (A2-I-A3) gZ, Ki = exp(fci), and 
iCi G KnA. Thus, we obtain {Adg{K)){a' -\-k") + 2ki = 
{Ads{K)){a' + k" + 2k2) = {Adg{K)){a' + ka) = tW, 
where Ki — exp(fci), Ki G Kd A and i G {1,2,3}. Corol- 
lary ^ completes the proof. □ 

Searching for a refinement of Cor. |21 we state bounds 
on the coefficients of ai, 02, and 03 of an element 
aiXr -\- a2Xs -I- 03X9 of a. It follows from Lemma|21and 
Lemma 13 that the coefficients a^, i G {1,2,3}, are peri- 
odic with period tt. (Concerning this periodicity, we refer 
also to Ref. 47, Appendix B] and Ref. |23, p. 7].) Bearing 
the TT-periodicity in mind, we can restrict the coefficients 
Oi to the interval [—^,^]. This choice is compatible with 
our conventions in Section IlII CI To reduce the symme- 
try induced by the Weyl group, we restrict ourselves to 
elements of the closed fundamental Weyl chamber. From 
Eq. 1)10(1 or from the s-order of Section HV Bl we get that 
fli > 0, 02 > 0, fli > 02, and 02 > a3. These considera- 
tions lead to the following corollary. 

Corollary 3 (see Ref. [6, Thm. 1] or Ref. pT, Re- 
sult 2]). Assume that H is a non-local Hamilton oper- 
ator acting on a two-qubit system and that we intend to 
simulate the unitary gate U. Let a and a' be elements 
of a, where a = aiXr -\- a2Xs -\- aaXg — {Adg{Ki)){H), 
a' = a[Xj + fljATg -I- 03X9, and U = K2 e-Kp{a')K3 for 
some local unitary operations Ki, K2, K3 G K. In ad- 
dition, we force ai, 02, as, a[, 02, and a'^ to be ele- 
ments from the interval [— f ,§]• We use the notation 
a = (ai, 02, 03)"^ and a' — {a[, 03, 0,3)'^ ■ 

A two-qubit system with Hamilton operator H and lo- 
cal unitary operators available is able to simulate the uni- 
tary gate U in time t if and only if the following equa- 
tion holds for at least one choice of z = (zi, Z2, z^)'^ G 
{(0,0, Of, (-1,0, OH; 

a' -\- TTz td. 

Remark. In the proof we follow Refs. 0,0- 

Proof. Due to Cor.|21it is sufficient to proof that for every 
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z € one of the following conditions holds: 



basis for two-qubit pure states: 



a +7r(0, 0,0)^ -<s a +7rz, 
a + 7r(-l, 0, 0)^ <s a + nz. 

We first consider the case that \zi\ > 1, for some i € 
{1,2,3}. Since a- < 7r/2, the maximal component {a' + 

7r£)i° of the s-ordered version of a' + ttz is greater than 
or equal to 2tt — tt/2 = 3tt/2. We check the conditions of 
Def . El and obtain that a' + 7r(0,0,Q)^ a' + ttz. 

Secondly, we consider the case that \zi\ < 1 for all 
i e {1,2,3}. By easy, but tedious, computations one can 
check that a' + 7r(0, 0, 0)^ <s a' + ttz for 

ze {(-1,-1,0)^, (-1,0, -If ,(0,-1,-1)^, 
(0,-l,lf,(0,0,Of, (-1,0,1)^} 

and that a' + 7r(— 1, 0, 0)"^ -<s a' + vrz for 

z e {(-1, -1, -If, (-1, -1, 1)^, (-1, 0, 0)^, 
(o,-i,of ,(o,o,-i)^,(o,o,in. 

For all other z £ {—1,0,1}^ we have that both a' + 
7r(0, 0,0)^ <s a' + TTZ and a' + 7r(-l,0, 0)^ a' + ttz 
hold. □ 



VI. LOWER BOUNDS FOR n-QUBIT SYSTEMS 

In the two-qubit case we used a particular decomposi- 
tion g = J + p of the Lie algebra which leads to a decom- 
position G = KAK of the Lie group where K = exp(t) 
is the set of local unitary operations. By this approach, 
e.g., the optimal simulation result of Fact HI can be ob- 
tained. In the more general n-qubit case we can use de- 
compositions g = t + p oi the corresponding Lie group 
where K = exp(t) contains all local unitary operations. 
Although K is in general not equal to the set of local 
unitary operations we can generalize the approach from 
the two-qubit case in order to prove lower bounds on the 
time complexity for gate simulation. Lower bounds were 
considered in Ref. |lOl. and we refine and generalize the 
approach of Ref. 10] in this section. In doing so, we put 
this approach in a broader context. 



A. Magic basis (for two qubits) 

We begin by recalling the Bell basis and the magic 
basis. The Bell basis (see Refs. ^HE3) is a vector space 



I1>- 



71 ^'^'^ 
71 ^'^^^ 



|10)): 
|10)), 



We employ the ket-vector notation, see, e.g., Ref. 0. 
If we include some relative phases in the Bell basis we 
get the magic basis which was introduced in Ref. |S| and 
coined by Hill and Wootters Q: 



ei) 






:=i|$- 




:= 


64) 





The magic basis is connected to the entanglement of for- 
mation, see Ref. jSj and related work in Refs. 50, 51]. We 
neglect here this connection, but refer to Section IVlIII 

The magic basis has two important properties. First, 
the local unitary operations on two qubits are real and 
orthogonal in the magic basis, see Ref. p. 5023] and 
Thm. 1 of Refs. j52]. Secondly, the elements of the A = 
exp(o) (for notations see, e.g., Sec. IIII 0|l are diagonal 
in the magic basis, as remarked in Ref. W^, p. 3] and 
Ref. 0, p. 2]. The basis change from the standard basis 
{JOO), JOI), JlO), 111)} to the magic basis is given by 
where 

/I i 
J_ I i 1 

V2 i -1 
\1 -i 

For elements U € SU(4) the map U ^ Q'^UQ (see 
Ref. HI) reflects the isomorphism between SU(2)(g)SU(2) 
and SO (4), see, e.g., Ref. ^ p. 52]. 



B. Representation theory 

It is not obvious how the magic basis generalizes to 
higher number of qubits and which properties remain. 
Motivated by the properties of the magic basis for two 
qubits, we seek for basis changes of the local unitary oper- 
ations (SU(2))®" into the orthogonal group (if possible). 
To analyze this we need some representation theory. 

Defin itio n 11 (Lie group representation, see, e.g., 
Ref. | 53L p. 210]). A complex representation of the 
Lie group G in the finite-dimensional and complex vector 
space Vc is a continuous homomorphism t: G GL(Vc) 
from the group G into the group GL(Vc) of invertible and 
linear transformations which operate on Vc. 
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A representation r in a finite-dimensional and com- 
plex vector space Vc is called irreducible if there exists 
no subspace Uc other than Uc = or Uc = Vc such 
that the subspace is T(G)-invariant, i.e., the equation 
r(G)Uc C Uc holds (see, e.g., Ref. 0, p. 210]). We 
state an important fact on tensor products of irreducible 
representations. 

Fact 9 ([54, Prop. 4.14 of Chap. II]). If n is an 

irreducible complex representation of Gi in the complex 
vector space Vc and T2 is an irreducible complex rep- 
resentation of G2 in the complex vector space Wc, then 
Ti®T2 is an irreducible complex representation of Gi x G2 
in the complex vector space Vc ® Wc . Furthermore, any 
irreducible representation of Gi x G2 is a tensor product 
of this form. 

Below we use bilinear forms B: Vc x Vc — > C, which are 
C-linear in both arguments, to characterize irreducible 
complex representations. Let vi and V2 be some arbi- 
trary elements of Vc. A bilinear form is called sym- 
metric if B{vi,V2) — B{v2,vi) and skew-symmetric if 
B{vi,V2) — ~B{v2,vi). A bilinear form is r(G)-invariant 
ifB{vi,V2) = B{t{G)vi,t{G)v2) for aU G e G. 



Definition 12 (cf. Refs. |53l l54|). Consider an irre- 
ducible complex representation t of G in Vc. The repre- 
sentation r is said to be of 

• real type if Vc admits a bilinear form which is 
nonzero, non-degenerate, T(G)-invariant, and sym- 
metric, 

• complex type if Vc admits no bilinear form which 
is nonzero, non-degenerate, and r(G)-invariant, 

• quaternionic type if Vc admits a bilinear form 
which is nonzero, non-degenerate, r(G)-invariant, 
and skew-symmetric. 

We introduce the map Xr: G ^ C,G 1-^ Tr(T(G)) 
which is the character Xt of the representation r. We 
use the character to characterize the type (real, complex, 
or quaternionic) of irreducible complex representations. 

Fact 10 ([53, Thm. 4.8.1]). Letr denote an irreducible 
complex representation of G in Vc . The character Xt is 
real-valued if and only if there exists a t{G) -invariant, 
nonzero, complex bilinear form B on Vc which is auto- 
matically non- degenerate and uniquely determined, up to 
a nonzero scalar factor. This bilinear form B is either 
symmetric or skew- symmetric. 

By means of Fact ^1 we can decide if the type of a 
representation is complex. To complete the classifica- 
tion of the type (real, complex, or quaternionic) of ir- 
reducible complex representations, we state another fact 
which allows to determine the type of a representation 
by computing an normalized integral over the compact 
Lie group G. 



Fact 11 (see Ref. fsi, Prop. 4.8.7] and Ref. [H 
Prop. 6.8 of Chap. II]). Let t be an irreducible com- 
plex representation of the compact Lie group G in Vc 
with character Xt ■ Then we have 



T is of real type, 
<^ T is of complex type, 
T is of quaternionic type. 



Representations can be identified with subgroups of 
GL(Vc), so we can extend Def. El to subgroups of 
GL(Vc). We denote the general linear group on a com- 
plex vector space of dimension k by GL(A;,C). Next we 
characterize the subgroups of GL(A:, C) that are conju- 
gated to subgroups of the orthogonal group (motivated 
by Sec. IVI A|) or the symplectic group. 



Fact 12 (adapted from Ref. |26|, Thm. H of 
Chap. 3]). A compact subgroup of the general linear 
group GL(A:,C) is conjugated in GL(fc,C) to a subgroup 
of the orthogonal group 0{k) if and only if it is of real 
type. Accordingly, a subgroup o/GL(2fc, C) is conjugated 
in GL(2fc,C) to a subgroup of the (unitary) symplectic 
group Sp(/c) if and only if it is of quaternionic type. 

Remark. Actually, Ref. [13] gives an algorithm to com- 
pute the basis change from the bilinear form mentioned 
in Def. ^1 For the notation Sp(A:) see Subsection IVI Dl 
and Ref. (SSj . 

After this preparation we consider the case of local 
unitary operations (SU(2))®". We employ the standard 
representation of SU(2): 



G 



a + ib c - 
— c -\~ id a 



id 
ib 



where a, 6, c, e M and + 6^ + + rf^ = 1. To compute 
the integral of Fact ^2 wc introduce the real parameters 
< (f> < 2tt, < tpi < tt, and < -02 < tt as follows: 

a = cos(0) sin(V'i) sin(V'2), 
b — sin(0) sin(7/;i) sm{ip2), 
c = cos{iJji) sin{ip2), 
d — cos{'tp2)- 



We obtain 




X.(G2)d(SU(2))«" 



Xr(G2)rfSU(2)j 



^((/),-!/'l,'!/'2)'iV'2dV'l# 

\0 

= (-1)", 

where 

S(0,'(/'i,V'2) =[4cos(</.)^(l - cos(V'i)^ - C0S(V'2)^ 

-I- cos('(/'i)^ cos(V'2)^) - 2] sin(^i) sin(-02)^- 
This proves the following theorem. 
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Theorem 3. The local unitary operations on an even 
number of qubits are conjugated to a subgroup of an or- 
thogonal group. The local unitary operations on an odd 
number of qubits are conjugated to a subgroup of a ( uni- 
tary) symplectic group. 

Remark. Similar results as in this subsection are obtained 
in Ref. jH^ using a different approach. 



C. Thompson's theorem and majorization 

Following Ref. i we present in this subsection a the- 
orem due to Thompson and a majorization condition 
for the spectra of the sum of two Hermitian matrices. 
Both results will be employed below. 

Fact 13 ([57]). Let A and B be Hermitian matrices. 
Then there exist unitary matrices Ui and U2 such that 



where x = (ii, . . . ,X2k) e C^*" and y = (yi, 



exp(iA) exp(i_B) 



exp(i[/f MC/i 



iUo^BU2) 



This result of Thompson relies partly on a conjecture 
of Hor n l58| . This conjecture was recently proven [sol 
i3, EUHl 113. By induction, we get the following 
corollary. 

Corollary 4. Let Aj denote Hermitian matrices. Then 
there exist unitary matrices Uj such that 



Yl exp(U,) = exp ^^J^J 



We state now a result which gives us bounds for 
the the spectra of the sum of two Hermitian matri- 
ces. Reference attributes this result to Ky Fan 
|65| . We denote the vector of eigenvalues of the k x 
fc-dimensional matrix A, including multiplicities, by 
spec(A) = (spec(yl)i, . . . , spec(A)fc)^. In addition, we 
assume that spec(A)i > spec(A)j if i < j (1 < i, j < k). 



Fact 14 (gll Thm. 9.G.I.]). Let A and B denote 
Hermitian matrices. Then the following equation holds: 

spec(A B) < spec(A) -|- spec(i?). 



D. Lower bounds 

In this subsection we derive lower bounds on the mini- 
mal time to simulate unitary operations (see Dcf ^ . We 
begin by discussing the (unitary) symplectic group. Fol- 
lowing Ref. p. 22], we introduce the bilinear form 



Bspix,y) 



E 

j=i 



i-\-k ^i-\-kyi) 



(16) 



-'2k 



, y2k) 



Let Jfc denote the matrix 



Ofe Ik 
-h Ofc 



where Ik is the k x fc-dimensional identity matrix and Ofc 
the k X fc-dimensional zero matrix. 



Definition 13 (see, e.g., Ref. | 55L Prop. 1 on 
p. 22]). The subgroup of the unitary group U(2fc) of 
degree 2k composed of the matrices M which leave the 
bilinear form in Eq. (|15|l invariant, i.e., which satisfy the 
condition 



JkM = Jfc, 



(17) 



is called the (unitary) symplectic group and is denoted 
by Sp(fc). 

The group Sp(fc) can be considered as operating on 
a fc-dimensional module over the quaternions H leaving 
a symplectic (scalar) product invariant |55i, pp. 16-24]. 
All elements of Sp(fc) have determinant one, see, e.g., 
Ref. |U p. 203]. When we regard Sp(fc) as a manifold 
its real dimension is 2fc^ -I- k ,5^ p. 23]. 

We recall from Thm. |31 that the local unitary opera- 
tions on an odd number of qubits are conjugated to a 
subgroup of a (unitary) symplectic group. Using that 
(Jk)^^ — —Jk, the condition in Eq. ((T7|l can be proved 
to be equivalent to 



M- 



JkM^iJk) 



(18) 



We know that the local unitary operations on an odd 
number of qubits meet the condition in Eq. 1181) in some 
appropriate chosen basis. But we can state the condition 
also in the standard representation of SU(2)'*" with 
odd. We use the identification 2fc 
the matrix 



2". Let J' denote 



1 

-1 



{iayf 



and recall the standard representation of SU(2): 

^ f a-\- ib c-\- id^ 
I — c -\- id a — ib J 

where a,b,c,d G M. We use the notation 

J \-Cj + idj Oj - ibj) ' 

where Oj, bj, Cj, dj G M. It can be checked that 
and we obtain that 



(19) 
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which holds obviously for n odd and even. From now on, 
n is no longer restricted to be odd. We emphasize that 
= {J'jr = (-l)"j;. It follows that 



(g)G, | j;, |(g)G, 



= j:. 



(20) 



Let Ti denote the 2"-dimensional complex vector space 
on which the group SU(2") operates. We introduce the 
bilinear form Bt-c{x, y) := JnV on the Hilbert space TL. 
We have that 

Bn{y,x) = y^J'^x = (-l)"x^42/ = {-irBn{x,y), 

which proves that Bh{x, y) is symmetric for n even and 
skew-symmetric for n odd. From Eq. (|20|l we get that 
Bn{x,y) is left invariant by SU(2)®". Hence, we have 
identified B-H{x,y) as the bilinear form of Def. 1121 oper- 
ating on Vc —li.. 

This motivates the following definition of the tilde 
mapping, which operates on the local unitary operations 
as the inverse operation. 

Definition 14. We introduce the tilde mapping 



SU(2") ^ SU(2") 



It is apparent from [JnU'^ {J'nY^J'nU'^ {J'n)~^]'^ = -^2- 
and det[( J4)~^C/"^ J^] = det(?7) that the tilde mapping 
preserves the group SU(2"). 

Remark. The tilde mapping is a generalization of the map 
U ^ {ay)®''U'^{ay)®'^ for even n from Ref. [H p. 5]. It 
can be easily checked that the two maps coincide in the 
case of even n. See also the discussion in Sec. IVIIII 

We state now an important lemma characterizing the 
tilde mapping. 

Lemma 4. Let V and W denote some local unitary oper- 
ations and let U denote some arbitrary unitary operation. 
The following equations hold: 

1. V = 

2. W = W-^, 

3. vuw-^ivuw) = vuuv-^. 

Proof. The first and second claim follows from Eq. (|19|l . 
We prove now the third claim: VUW'^iVUW) — 
VUWWUV ^VUWW-^UV-^ = VUUV''\ □ 



Theorem 4. Assume that H is a non-local Hamilton 
operator acting on an n-qubit system and that we intend 
to simulate the unitary gate U . An n-qubit system with 
Hamilton operator H and local unitary operators avail- 
able is able to simulate the unitary gate U in time t only 
if the following equation holds for at least one choice of 
ze 1?" : 



arg 



spec(t/J7) + 27rz 2tspec(ff) 



Remark. This theorem generalizes the work in Ref. [T3, 
Thm. 5 and Cor. 7]. In the proof we use ideas from 
Ref. [13. 

Proof Assume that U = [Hjli exp(-itji7)] Wq is a 
simulation of U, where Wj denotes some local unitary 
operations, tj > 0, and = ^- Let 



y ^ I'^rn for j = m, 
' \Vj+iWj forO<j<m. 



and Hj = VjHV^^ (1 < j < "m). Thus, the simulation 
of U can written as [/ = [Iljli 6^P(~*^j-^j)]^ where 
spec(i?j) = spec(iJ) for all 1 < j < m. 

We use Cor. 0] to find some suitable Hermitian oper- 
ators -ffj and H'- with spec(iJj) — spec(iJ") ~ spec(_ff) 
and j e {1, . . . , m} such that 



J|exp(-itji?j) 

— exp 



[]cxpHt,iJ,) 



(21) 



When we combine Vq — Vq ^ with Eq. (|21() we obtain 
that 



spec(C/C/) = spec < exp 



iY.t,{H'^+H^) 



We employ Fact to complete the proof: 



arg 



spec(C/[/) 



+ 2'Kz = spec 



-< 2tspec{H). 



□ 



This proves that local unitary operations preserve the 
spectrum of UU. We state now the theorem which gives 
us lower bounds for the minimal time to simulate a uni- 
tary gate. We use the notation arg where arg [exp(irn)] = 
m and arg [(xi, . . . ,xi)'^] = (arg[a;i], . . . , arg[a;i])^. 



E. Involutive automorphisms 

We end this section by highlighting connections be- 
tween the tilde mapping of Def. 1141 and involutive auto- 
morphisms of the Lie algebra su(2"). 
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The tilde mapping is similar to the ()*-map used in 
the proof of Lemma |21 in Section ]V\ For n odd K must 
be equivalent to Sp(2"^^) and respectively for n even K 
must be equivalent to S0(2"). In both cases, the map 
U UU plays a similar role as the map (j> in Section Ivl 

Following Ref. [2^ pp . 451-452] we state all Rieman- 
nian symmetric spaces SU(2")/i^' induced by involutive 
automorphisms of the Lie algebra su(2"). We have to 
consider three cases which correspond to the types AI, 
All, and AIII of involutive automorphisms. In the case 
of type AI we have to treat the Lie algebra g = su(fc) 
and the involutive automorphism 9ja{g) = g* . The in- 
volutive automorphism O^i gives rise to the Riemannian 
symmetric space SU(fc)/SO(fc). 

The Lie algebra g = su(2fc) and the involutive auto- 
morphism 6'aii(5) — Jkg*{ Jk)~^ belong to type AIL We 
obtain the Riemannian symmetric space SU(2fc)/Sp(fc). 

For completeness we mention the type AIII even 
though we do not use the corresponding Riemannian 
symmetric space in this paper. The Lie algebra is 
g = 5u{p -\- q) and the corresponding involutive auto- 
morphism is given by 6'aii(5) — Ip,q9lp,q- We have used 
the notation 

T _ ( Op.? 1 

where Ip denotes the p x p-dimensional identity ma- 
trix and Op^q denotes the p x g-dimensional zero ma- 
trix. This gives us the Riemannian symmetric space 
SU(p -t- g)/S(U(p) X U(g)). The group S(U(p) x U(g)) 
can be represented by the matrices 

f 91 Op^q\ 
\%,P 92)' 

where gi G U(p), 52 e \J{q), and det(5i) det(g2) = 1- 

VII. RELATED WORK 

To recognize the considerable amount of related work 
we give a short outline of the connections to our work. 
Various aspects of (infinitesimal) Hamiltonian simulation 
as considered in Sect . II V I we re studied in Refs. ^ [isL IT^ 

iiiiiiiiiiiiiiiiaiiiaiMiMizaEipizliM 

Some references consider models where — in contrast to 
our model — additional resources were used: prior en- 
tanglement [T^I. additional classical communication , 
measurements ]^ : ancillas 0, . 

For two qubits, gate simulation (see Sect. 0) was an- 
alyzed in Refs. H II S 111 El El- In Ref. Lie 
group decompositions were used to obtain a theory of 
n-qubit gate simulation. In general these decomposi- 
tions do not lead to optimal simulations. In the case of 
three qubits some progress on the time optimality prob- 
lem f or g ate simulation was reported in Ref. |77| , see also 
Ref. |78j . Concerning lower bounds, we have generalized 
(see Sect. IVlj) the approach of Ref. US ]. 



VIII. DISCUSSION 

In this section we address a peculiar similarity between 
our approach to lower bounds on the time complexity for 
gate simulation and the concurrence Ifl, [III ISU , as 
well as some of its generalizations [s^ IsillS^ l83l l84 Issl 
m HI m 189', 90, 91 m 111 Hi . The concurrence C of 
a pure two-qubit state & C'* was defined in Ref. 
as 

= 1(^1^)1, 

where :— {ay (g) ay)(\'ip))*. Let Ai, A2, A3, and A4 
denote the (positive) square roots of the eigenvalues of 
the matrix p/j, where p :— {ay ® ay)p*{uy ® OjJ. We 
assume that Ai > A2 > A3 > A4. References H0 show 
that the concurrence C of a two-qubit density matrix p 
is given by 

C{p) = max{0, Ai - A2 - A3 - A4}. (22) 

Uhlmann [s^l considered generalizations of the con- 
currence. Following this approach we introduce some 
notations. Let us call a map ^ that operates on 
a complex vector space Vc antilinear if the equation 
i?(6i|V'i) + b2\^2)) = {biy'&{\iPi)) + (62)*i?(|V'2)) holds 
for all 61,62 G C and all jV'i), |V-'2) G Vc- For an anti- 
linear operator d the (Hermitian) adjoint i?^ is defined 
by the condition that (?Ai|i5^'02) — {4'2\'&'4'i) holds for 
all \tpi),\4'2) G Vc. If an antilinear operator satisfies 
the condition d'^ = 'd^^ we call this operator antiuni- 
tary. When the map d is antiunitary and = d holds, 
then we have that 1?^ equals the identity map and we 
define ■(? to be a conjugation. A skew conjugation is an 
antiunitary operator d fulfilling 'd~^ = —d. Assume in 
the following that is a conjugation. Now, Uhlmann 
defined a generalized tilde mapping by its operation on 
pure states It/i) := •&{\i')) and its operation on density 
matrices p :— dpd^^ = dpd. In addition he generalizes 
the concept of concurrence to more than two qubits for 
pure states 

c^{m 

and for mixed states 

C^{p) := min^|(0j|^j)|, 

j 

where the minimum is taken over all decompositions p = 
X]j \4'j){4'i\ of P into non-normalized pure states \4>j)- 
Uhlmann [s^ proved that in strong analogy to Eq. I|22l) 
the generalized concurrence is given by 

C{>{p) = max < 0, Ai - Aj > , 
I J 

where the A^s are the square roots of the eigenvalues of 
the matrix pp and Ai > Xj for i < j. 
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References [s^ |^ |^ consider the map on density 
matrices given by 

p^p:=K)«VK)^"- (23) 
In addition, the map 

iH ^ iH (-iaj;)®"(ii7)* {{-icTyf")'^ (24) 

is introduced in Ref. for elements iH of the Lie al- 
gebra 5u(2"). The map in Eq. (Pljl can be apphed to a 
Hamilton operator H: 

H = i-iayr-iHy ii^iayr-)-' = i^ryr^iHY {ay)^- . 

This shows that both Eq. 123) and Eq. (|23I) are induced 
by the conjugation §1 given by 

In this case we get the concurrence C^}^. The corre- 
sponding tilde mapping is given by its action on pure 
states = and its action on density matrices 

p = §ipdi~^ = dip^i. A result of Ref. Prop. 8] (for 
related remarks see Ref. jS^i) states that the conjugation 
•di is the (up to a phase) unique antilinear operator act- 
ing on the complex vector space (C^)®" which is invariant 
under basis changes by local unitary operations U except 
for an factor equal to det(?7). Further on, Ref. ^| states 
that such an antilinear mapping exists only for n-qubit 
systems, and not for general n-qudit systems. 

After this short excursion into entanglement measures 
of concurrence-type we can state a connection between 
this type of entanglement measures and lower bounds 
on the time complexity for gate simulation. The tilde 
mapping of Def. 1141 which was used in the main body of 
the text, can be interpreted in the context of concurrence- 
type entanglement measures. Since i/^ = H* holds for 
all Hermitian operators we obtain 

*([/) = j;(exp(ii/))^(4)-i = exp(i4i/*(j;)-i), 

where iJ'^H* (J^^)^^ is up to a minus sign equal to the rhs 
of Eq. (|24|l . This highlights that if we consider the lower 
bounds introduced in Sec. IVII we are essential in setting 
of Ref. m with 1? = 

In both constructions, for the lower bounds on the 
time complexity to simulate unitary operators and for 
the computation of the concurrence, the essential point 
is that the spectrum of both UU and pp is invariant un- 
der local unitary operations. This highlights that there is 
a connection between entanglement measures and lower 
bounds on the time complexity for gate simulation. We 
are looking forward to generalizing some of these ideas. 

IX. CONCLUSION 

Starting with an extensive reconsideration of infinites- 
imal Hamiltonian simulation and gate simulation in the 



two-qubit case, we streamlined the different approaches 
by using Lie-theoretic methods. As the success of this 
approach suggests, this seems to be the appropriate level 
of description for such a theory. 

Going beyond two-qubits, we derived lower bounds on 
the time complexity for gate simulation. For this aim, 
we developed an analogon of the magic basis for general 
multipartite qubit-systems. This gives us a first idea of 
the structure of unitary operations w.r.t. the set of local 
unitary operations. In addition, we related our approach 
to entanglement measures of concurrence-type. 
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APPENDIX: SPECTRAL THEORY FOR 
INFINITESIMAL HAMILTONIAN SIMULATION 

In this appendix we prove a similar version of Thm. |21 
Relying on Subsection IIV Bl we use arguments from the 
theory of majorization applied to the spectrum of Hamil- 
ton operators. We denote the vector of eigenvalues of the 
k X fc-dimensional matrix A, including multiplicities, by 
spec(A) = (spec(A)i, . . . , spec(A)fc)^. In addition, we 
assume that spec{A)i > spec(A)j- if i < j (1 < i,j < k). 
The majorization of the spectra of two matrices is, by a 
theorem of Uhlmann related to the convex combi- 
nation of unitary orbits. 

Fact 15 (Uhlmann, see, e.g., (9^ Satz 3] or [i^, 
Thm. 2-2.]). For Hermitian matrices A and B the con- 
dition spec(j4) ~< spec(i?) is equivalent to 

A = J2q^Ur^BU,, 

i 

where Ui are unitary matrices, qi > 0, and X^i'Zi = 1- 

We need another theorem connecting the notion of ma- 
jorization with the convex hull of all permuted versions 
of a vector. 

Fact 16 (Rado, see, e.g., Ref. \9T) or Ref. fill 
Prop. 4.C.I.]). The vector x is majorized by the vector 
y if and only if x lies in the convex hull of all permuta- 
tions of y. 

The spectral version of Thm. |2] reads: 

Theorem 5 ([46, pp. 9—10]). Assume that H and H' 
are non-local Hamilton operators acting on a two-qubit 
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system. Let a and a' be elements of a, where a = aiX^ + 
02^8 + ogXg = (Adg (L))(iJ) and a' = a^X^ + a^Xg + 
a'^Xg = {Adg{L')){H') for some local unitary operators 
L and L' . 

A two-qubit system with Hamilton operator H and local 
unitary operators available is able to simulate the Hamil- 
ton operator H' in time t if and only if 

spec(a') ^ t spec(a). 

Remark. The necessity of this condition was proven in 
Ref. [13. In the prove we foUow Ref. 

Proof Due to Fact □ we can choose a and a' as given. 
The "only if" -case foUows by Fact El We consider now 
the "if" -case. By invoking Fact 1161 we get that 

spec(a7i) =^qkPkSpec{a), 

k 

where Pk is a permutation, qk > 0, and qk — 1. Since 
a and a' are elements of a, they commute. It follows that 
there exists a basis where a and a' are simultaneously 
diagonal. In that basis the permutations Pk correspond 



to permutations of the diagonal elements of a. For that 
reason we have that 

{a'/t)=Y,QkU;:^aUk 

k 

for some unitary operators Uk which permute the spec- 
trum of a. We emphasize that the Uk are not neces- 
sarily local. But we prove now that we can find lo- 
cal unitary operators implementing any permutation of 
the spectrum of a. Conjugation by the local unitary 
operators ((ctq - icri)/\/2) (g) ((ctq - *o'i)/V2), ((cro + 
ia3)/V2)(E)iiao+icr3)/V2), and (((Tq -hiai)/\/2) «) ((ctq - 
iai)/V2) permutes the eigenvalues as respectively fol- 
lows: (1,2,3,4) ^ (2,1,3,4), (1,2,3,4) ^ (1,3,2,4), 
and (1, 2, 3, 4) 1-^ (1, 2, 4, 3). As all permutations on four- 
vectors are generated by this permutations, the "if" -case 
follows. □ 

We note that the local unitary operators that permute 
the spectrum of elements of a are given in Ref. |4^ , but 
there the second local unitary operator is misprinted. 
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